Neutron reflection from condensed matter, the 
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Abstract 

The Goos-Hanchen (G-H) effect for neutron reflection from condensed matter is 
considered. An experiment to quantify the effect is proposed. The relation of G-H 
shift to the neutron coherence length is considered. 
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packets. 
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1 Introduction 



Newton suggested first (see, for example [1,2]), 
that a beam of light at a specular reflection shifts 
some distance along the surface of a reflecting 
mirror, as shown in fig. 1. If the incident light is 
represented by a ray, then propagation of the ray 
in a mirror can be described by some trajectory, 
where exit point B of the reflected ray does not 
coincide with entrance point A of the incident 
one. The shift AB is called G-H effect, because F. 
Goos and H. Hanchen [2,3] experimentally mea- 
sured it. 

The shift should exist not only for light but also 
for particles, because propagation of particles in 
quantum mechanics is described by wave func- 
tions similar to the wave field of the light. 

The shift depends on interaction of particles with 
matter and investigation of it can be useful for 
applied [1,4] and fundamental research. In this 
article we consider G-H effect in neutron physics 
(see, for example, [5] and references there in). 
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We address several problems. First, we calculate 
longitudinal G-H shift at total reflection for re- 
stricted Gaussian beams (like in optics) and for 
wave packets. We believe that this approach is 
more general than the one (with the help of two 
plane waves) used in [5]. Second, we calculate de- 
viation of reflected beam from specular direction 
for Gaussian beam (it is similar to [6] in x-ray op- 
tics) and for wave-packets. And third, we discuss 
the problem of enhanced G-H effect at reflection 
from a thin layer evaporated on a totally reflect- 
ing substrate. We believe that investigation of 
G-H effect in optics of neutrons can help to un- 
derstand deeper the nature of the neutron wave 
function: whether it is a plane wave or a wave 
packet. If it is the wave packet then whether it is 
related to preparation of the neutron beam or it 
is an intrinsic property of the neutron. 

Section 2 reminds the reader how the G-H ef- 
fect is calculated. We consider plane waves, finite 
beams and wave packets reflection from a single 
interface and from a layer on a substrate. In sec- 
tion 3 we discuss how peculiarities related to the 
G-H effect can be measured, and what informa- 
tion on neutron wave packet can be obtained. We 
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Figure 1. Schematical explanation of the G-H shift 

address there a question whether calculation of 
reflection amplitude from a layer with two inter- 
faces, which uses coherent superposition of co- 
herent multiple reflection of waves (see, for ex- 
ample [7,8,9,10]), is accurate, and justify the ap- 
propriateness of introduction in [11] of so called 
"coherence length" . This length determines such 
a length of the G-H shift, at which reflection co- 
efficient contains a considerable contribution of 
incoherent reflection. 



the incident particle. Below we consider three 
types of this function: plane wave, stationary 
beam, and a wave packet. 



2.1 Plane wave 

Let the incident particle be described by a plane 
wave 

exp(ifco^" — icut), 

where ct; = A;q/2 is the energy of the neutron and 
fco = —kz) is a wave vector, with component 
parallel, and — normal to the interface. 

The reflected wave is 

Rexp{ikzZ + ik\\r\\ — iut) = 
\R\ exp{—2i(f){kz) + ik^z + — icot), (1) 



At the end of this paper, in conclusion, we shortly 
discuss applicability of our calculation of the neu- 
tron G-H effect to x-rays 



where the reflection amplitude 



R{k,) 



kz k'^ 
kz + k' 



2 Calculation of the G-H shift 



{k'^ = s/k'^ — u) is represented as i? = \R\ exp(— 2i0). 



Here we consider neutron reflection from an ideal 
surface of a semiinflnite mirror, which means that 
we have a single interface. Reflection takes place 
because interaction of neutron with matter. This 
interaction is characterized by an optical poten- 
tial m/2, where u = ATiNob, Nq is atomic den- 
sity, b is coherent scattering amplitude on a sin- 
gle atom, and for convenience we set h = m = 1. 
If normal component kz of the incident neutron 
momentum k is small {k'^ < u), the neutron is 
totally reflected from the interface, i.e. reflection 
coefficient \R\'^ = 1. If normal component kz is 
large {k^ > u) , the reflection is not total (in the 



following we call it "nontotal" ) , and 
Here R denote reflection amphtude. 



\R\' < 1. 



Reflection coefficient depends not only on k and 
u, but also on structure of the wave function of 



The phase (j) and the absolute value \R\ depend 
only on kz- When k^ > u and m/2 is purely real, 
i.e. there are no losses, then \R\ < 1 (reflection 
is nontotal), and 4> — 0. Since there is no phase 
shift = 0, the phase of the reflected wave at 
every points on the interface is equal to the phase 
of the incident wave. 



When kj < u then \R\ = 1, i.e. reflection 
is total, and the phase is nonzero: < (f) — 

arccos(A;^/y^) < n/2. It means that the phase 
of the reflected wave (1) is equal to the phase 
kzZ + fc||r|| — ujt of the incident wave at points 
ry displaced along the interface by a distance 
6r\\ — 2(j)/k\\. This distance can be called the 
G-H shift, but a spatial displacement can be 
well deflned only for flnite size beams. In case of 
inflnite plane waves we can talk only about the 
phase shift. 
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2.2 Two plane waves 

In [5] the G-H shift was calculated with the help 
of two plane waves. The incident neutron was 
supposed to be described by the wave function 

■01 = exp(ifeir) + exp(ifc2r") 

, consisting of two waves with wave vectors fcj = 
{kix-i ki±) , where kx is a component along surface, 
and /cj_ is perpendicular to it. This function can 
be also represented as 



■01 = 2exp(ifcr) cos{Skr), 



(2) 



where fe = (fcs + fei)/2, 5k = (fci - k2)/2. We 
see, that at the interface ^ = we have a plane 
wave ex.p{ikxx) modulated by cos{Skxx). 

After total reflection this wave function is trans- 
formed into 

t/jr — R{ki±) exp(ifc^r) + R{k2±) exp(ifc2r) = 

= exp(ifcir + + exp(ifc2^ + i(f)2), 
where for total reflection we replaced reflection 
amplitudes R{ki±), which depend only on nor- 
mal komponents of the wave vectors by exp(i0j). 
Wave vectors of the reflected waves are kl — 
{kix, —ki±). The sum of the reflected waves can 
be also represented as 



ijjr — 2exp(ife''r -|- icj)) cos{5k^r + 54>), 



(3) 



where fc" = (fc^ + k[)/2, = (0i + 02)/2, 6k = 
{k[ - fc^)/2, and (50 = (01 - 02)/2. We see, that 
at the interface z — we have a plane wave 
exp{ikxX + (()) modulated by cos{SkxX + 5(l)). If we 
represent 6(f) as 6kx^ (why not?) we obtain that 
modulation of the reflected wave is shifted along 
X axis by ^ = 6(f)/ 6 kx, which can be considered 
and was accepted in [5] as the G-H shift. 



2.3 Wave function of a beam 



Imagine now that the wave function of a neutron 
is confined in space around a ray, i.e. it looks like 



a cylindrically symmetrical Gaussian [6] . Such a 
wave function can be stationary, and it creates an 
eUiptical spot at the interface [6] . We will place 
an origin r = of the coordinate system in the 
center of this spot and choose xz to be the plane 
of incidence. The beam at the interface is 



exp 



■0(r|| ,z — 0,t) — exp{ikxX — icut) x 



cos^ 6 y'^ 

+ — 

2 2 



(4) 



where r\\ = {x,y), uj = k'^/2 = {kl + kl)/2, 
1/s is the beam's width, and cos 6' = kz/k = 

^Jk"^ — k^./ k is the effective angle of incidence. 

Once we find the the wave function of the re- 
fiectcd beam at the interface, the position of the 
center of the reflected beam with respect to the 
origin will give the magnitude of the G-H shift. 

To find the refiected wave function we need first 
to restore the complete incident one from its spot 
(4). In order to do that we represent (4) as a 
Fourier expansion 



V'o(r-||,-2 = 0,i)= j 



27rs2 cos 9 



X 



exp(iq||r|| — iujt) exp 



(?a; kxY' 

2s^ cos^ 9 



,(5) 



and multiply the integrand by a factor exp {—iq^z). 
Then, 



exp 



{Qx kxY' Qy 



2s^ cos^ 9 s 



2 ' 



(6) 



where qz must satisfy q^^+ ql — 2uj — k"^ so that 
the function (6) satisfies the free Schrodinger 
equation. Thus = ^k'^ — q^y 

The incident wave function can be represented as 
a superposition of plane waves, while the refiec- 
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tion of a plane wave is known. Thus, we can im- 
mediately find the wave function of the reflected 
particle: 



^r{r\\,z,t) 



27rs^ cos 9 

oo 



2.3.1 Nontotal reflection 

At > u the cross-section of the reflected beam 
at the interface z — 0, according to (7), is 



V'r(r'||,^ = 0) = j 



27rs2 cos 9 



X 



exp 



{Qx kx)^ Qy 



2s^ cos^ 9 



(7) 



where R{qz)is the reflection amplitude of the in- 
cident plane wave exp(iqr||r|| — iqzz). This ampli- 
tude is well known: 



i?(g.) = ^^ = exp(-2x), 
Qz + q'z 



(8) 



exp 



exp(iq||r|| + {q^ - kx)Qx 



Is^ cos^ 9 



(11) 



where for simplicity we omitted the time factor 
exp(— ia;t). After integrating over d^gn we obtain 



V^^(r||,^ = 0) = i?(A;,)expf— C^^os^^j x 



where x =arcch(q'2;/-y/w) =arcsh(g^/y^), and 

<iz 



Since qz is not an independent variable, 



k"^ — q1, the amplitude R through depends 
orLq»:R[qz) = exp(—2x(g)), where g = |q|||,and 



x(g) = arcsh(\/ k"^ — q^ — u/\/u). 



(9) 



We assume for simplicity that s is small, i.e. the 
radius 1/s of the beam is large. In that case 
the integral (7) can be easily calculated asymp- 
totically. The function 2x(g) in (7) is expanded 
around maximum point q = k\\ = {kx, 0) of the 
Gaussian up to the linear term 2x(g) = 2xo — 
{q-k\\)^, where Xo =arcsh(A;'z/ , ^ = (^^, 0), 
and 



d 



= 2- 



(10) 



When k1 < u (total reflection) the wave-vector 
component k'^ and thus become imaginary, 

K = iK = - kl and = -ilx = 
likx / kzk^ . 



exp I - Y ^ + y^] \ X 



exp(i(A;a; + ^s^ cos^ 9)x). 



(12) 



One can see that there is no G-H shift, since the 
reflected beam is centered at the same point as 
the incident beam. However the reflection am- 
plitude is a little bit larger than for plane waves 
by the factor exp ( cos^ 9^ , and one can find 
that it is not specular. The wave vector compo- 
nent kx of the reflected beam is larger than that 
of the incident beam by the amount 



5kx = cos^ 9^2 



h.2y 



(13) 



and because of the energy conservation, the com- 
ponent kz must be smaller than that of the inci- 
dent wave. Thus the reflected beam shghtly turns 
from specular direction toward the interface. 

Deviation of the reflected beam from specular 
direction is easily understandable, if we take into 
account that the reflection coefficient is smaller 
for larger g^. Thus the refiected beam is enriched 
with smaller g^ (and larger g^ because of energy 
conservation). Thus the average q of the reflected 
beam deviates from the specular k. 
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2.3.2 Total reflection 

At k1 < u the section of the reflected beam by 
the interface according to (7) becomes 



V'r(r-||,^ = 0) = y" 



COS 9 



X 



exp 



2s^ cos^ 9 ^ 
exp(iq||r|| -i{q^ - kx)ix), 

where 



X 



(14) 



(15) 



Integrating over d^gy we obtain 

%l)r{r\\,z — — R{kz) exp{ikxx) x 

ew(-'^[{x-Lrcos'9 + y']\. 



(16) 



We see that the spot center of the reflected beam 
is shifted with respect to that of the incident one 
by a distance which is the G-H shift. The re- 
flection in this case is completely specular. How- 
ever, the phase k^x is identical to that of the inci- 
dent wave, i.e. it does not contain the correction 
^x^x^ which we would expect because of the G-H 
shift. 

Above, we limited ourselves to the linear term in 
the expansion of the exponent x iii (§)• If we re- 
tain quadratic terms we can find the broadening 
of the reflected beam. This, however, does not 
bring any new insight because of the Gaussian 
beam own broadening. 

Expansion of xiQ) is valid only when \k'^\ 3> s. It 
is no longer valid near k'^ = because C,x diverges 
at this point. If we want to find the reflected 
beam for the critical point k'^ = 0, we have to 
approximate exponent x^J ^ function 



l2kx{qx 



u 



(17) 



and calculate the integral (7) using the steepest 
descent method. As a result the spot of the re- 
flected beam at the interface is deformed, and the 
reflection is not exactly specular. This case will 
be analyzed elsewhere, when it will be clear hot 
to study properties of the single neutron wave 
function experimentally. 



2.4 The G-H shift for a wave packet 



In the previous paragraph we consider the G-H 
shift for a particle represented by a stationary 
bcam-likc wave function. It is logical, however, to 
describe a particle by a moving wave packet. The 
wave packet is not a stationary wave function, 
and its Fourier representation is 

7po{r, = y Mq) exp(igr - iujqt)d^q, (18) 



where all components of the vector q — (qpQz) 
are independent variables, and ujq = q'^/2. For 
coefficients A{q) we use Gaussian function. 



(27rs2)3/2 



exp(-(g-fc)V2s2), (19) 



where for incident particle k = (fcy, g^) 

Integrating (18) over q gives the Gaussian wave 
packet in space 



^o(r',t) 



(1 + isH)^/^ 



X 



exp 



2{1 + isH) 



{r - kotf + ikor - iukt , (20) 



where s'^t characterize spreading of the wave 
packet. We suppose that s is small and neglect 
spreading. 

Cross-section of the wave packet at the interface 
depends on time. We can choose the time when 
the wave packet center crosses the interface, t = 
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0. The cross-section at that time is 



^(r||,^ = 0,t = 0) =exp rj + ife||r||j .(21) 



The reflected wave function is 



i^rir, t) = j R{q,)A{q) exp{iq^r - iujqt)d^q, (22) 



where = (qpQz) {Qz is independent variable) 
and R{qz) is given by (8). 



because of motion of the reflected wave packet 
its center crosses the interface. It happens at the 
moment t — Cz/kz, and at this moment the cen- 
ter of the wave packet is located at r\\ = k\\(z/ kz- 
This location is shifted with respect to center of 
incident packet. So this shift well corresponds to 
the intuitive understanding illustrated in fig. 1, 
and can be called the G-H shift (15). 

It is important to note that the cross-section 
of the wave packet at the interface has a circu- 
lar form, not an elliptical form as in the case 
of a beam. Since the linear expansion is valid 
for k'^ > s, the G-H shift is ^ -C 2/s, and 
^11 = 2k\\(z/kz ^ 2k\\/skz respectively. 



2.4-1 Total reflection 
First we consider the case of the total reflection 

kl < u, when x in (8) is = arcsin( Jl - g2/„). 2.4.2 Nontotal reflection 



We assume that s is small and expand 2(j){qz 
near qz = kz up to the linear term, 2(f){qz) = 
2(f){kz) - Cz{qz - kz): where 



C. = 2/k': = 2/Ju - kl 



(23) 



Placing R{qz) = R{kz) e^{i{qz-kz)Cz) into (22) 
and taking into account (19), we obtain 



Mr^t) fti R{kz) j A{q) x 

exp(iqf^r + iC,z{qz - kz) - iujqt)d^q. (24) 
Integrating over q (neglecting spreading) gives 
^pr{r,t) ^R{kz)x 



expi-—(r + C-kt)^ + ikr-iu;kt\, (25) 



where k = (fcy, k^). and the vector ^ has compo- 
nents C = (0, 0, C). 

At the moment t — 0, when the incident wave 
packet crosses the interface, the center of the re- 
flected wave packet is at the point z = —(z under 
the surface (depth C in fig. 1). Thus, the real shift 
is not in x-direction. It is in 2;-direction. However 



We consider now the case k^ > u, when 

X —&xcsh.{^q^/u — 1) is real. Expanding 2x{qz) 
near = kz up to the linear term gives 
'^xi.Qz) = 2 x{kz) + Cziqz - kz), where = 
2/k'^ = 2/^k'^ — u. Then we place R{qz) = 
R{kz) cxp(— (^2 — kz)C,z) into (22), take into ac- 
count (19) and integrate the result over q, while 
neglecting spreading of the wave packet. As a 
result, the wave function is 



il^rir.t) = R{kz)e 



X 



exp 



where k' 
k^/2. 



(r — k'tY -\- ik'r — iujk't , 



(26) 



{k\\,kz - ), and ujk> = k'^/2 < 



We see that the center of the reflected wave 
packet crosses the interface at the same moment 
t = and at the same point ry = as the in- 
cident wave packet, i.e. there are no G-H shift. 
However, the momentum of the reflected wave 
packet is less than that of the incident wave 
packet. The normal component of the momen- 
tum is less by the amount 



Sk^, 



(27) 
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which means non-specular reflection. The energy 
cuk' of the reflected wave is less than that of the 
incident wave. However, since the total energy 
should conserve after the elastic interaction with 
the interface, the refracted wave has larger energy 
than the incident wave. 

The change in energy is the result of the indepen- 
dent reflections of the different plane waves, and 
thus the result of the deformation of the wave 
packet at the reflection. The higher is the nor- 
mal component of the component plane wave, the 
less is its reflection amphtude. To avoid such de- 
formation we must assume that the wave packet 
is an intimate property of the particle as, for 
example, is in the case of the de Broglie wave 
packet [12]. For the de Broghe wave packet the 
energy after the reflection does not change, and 
the decrease of the normal component of the 
wave vector means that the direction of the re- 
flected wave packet is different from specular one. 



2.5 G-H shift for reflection from a layer on a 
substrate 




Figure 2. Reflection from a layer of thickness d 
with the potential ui, placed on a substrate with 
the potential U2 > ui. For a sufficiently thick first 
layer the G-H shift AB can be arbitrary large. 

Calculations of the shift for a beam-hke wave 
function or for a wave packet for the mirror in flg. 
2 are the same as before. The only difference is in 
the form of the reflection amplitude in (7). It can 
be calculated with the help of the method [10], 
which takes into account multiple reflection of 
plane waves from the two interfaces. 



[1 - r^^)ri2exp{2iq^id) 
1 - riori2exp (2^5^10?) 



rpi + ri2exp(2ig^iri) 
1 - riori2 exjp{2iqzid) ' 



(29) 



The G-H shift for the total reflection is usually 
much smaller then the dimensions of the cross- 
section of the wave function at the interface. 
Since we assumed that k" s, then 



h h" 



«;2-^ 



k^s 



s cos 9 



(28) 



The shift becomes considerably larger, if the mir- 
ror consists of two parts as shown in flg. 2. It 
could be, for example, a layer of Al evaporated on 
Be substrate. We denote the potential of the top 
layer to be Ui and the bottom U2, and we choose 
Ml < U2. We also choose, that the total reflection 
takes place at the second interface between the 
layer and substrate. We can imagine, however, 
that the total reflection happens instead at the 
flrst interface of the mirror, only this reflection 
happens with a large G-H shift and with deep 
propagation (comparable to the thickness of the 
flrst layer) into the mirror. 



where 



qzi 



Qz + Qzi 



ri2 



Qzi - qz2 
Qzi + qz2 ' 



(30) 



qzi,2 = \jql-ui^2- 

We consider the case when q1 - 
at the same time U2 — q1 ^ s 
the reflection amplitudes from the substrate and 
from the whole mirror are unit complex num- 
bers (the losses are neglected). If we denote ri2 
exp(-2 i(/)i2), w here 0i2 
Q.Z2 = \lu2- ql, then 



- ui ^ s and 
^ In this case 



arcsin ( q'^^ / \Ju2 — ui) , 



R{qz) = 

ad-2i4„2 1 + ^01 exp(-2zg^id + 2i0i2) 



1 + roi exp(2iqzid - 2i0i2) 
exp{2iqzid - 2i0i2 - 2i0o2) (31) 
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where 0o2 



arcsin 



rpi sin(2g;;id - 2^12) 



^ Y^l + rgi + 2roi cos(25^id - 20i2) ^ 



.(32) 



In the case of a beam-hke wave function, expand- 
ing the phase 2qzid — 20i2 — 20o2 in (31) over 
qy —k\\, where k\\ = {k^, 0), up to the hnear term 
gives the G-H shift, 



+ 



'2k cr 



1 + rpi cos{2q^id - 2^12) 
1 + 2roi cos{2qzid - 2^12) + r^i ' 



X 



(33) 



The first term of the expression in the brackets is 
kzd times larger than the shift at a single inter- 
face. The factor outside the brackets depends on 
the neutron wave length and varies in the range 
from 1/(1 + roi) to 1/ (1 - Tqi). 

Same result can be obtained in the case of a wave- 
packet wave function. 



3 Possible experiments to measure the 
neutron G-H effect 



In this section we estimate the value of the G-H 
effect and look into perspectives of its measure- 
ment. 



3. 1 Nonspecular reflection from a single inter- 
face 



In case of a nontotal reflection of a beam, 
the component k^ — ksinO {6 is the angle 
of incidence) of the wave vector in the re- 

ficction plane along the interface increases by 
2kxkzs'^/k'^k'^{13). This is equivalent to rotating 
of the wave vector by a small angle 7: 



k sm.{6 + 7) = /Cx cos 7 + kz sin 7. 



From (34) follows that 5kx = kz'j, and when com- 
paring the latter with (13) we find 
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) 

'1-21./ ' 



(35) 



(34) 



In order to estimate 7 we need to find s. 

The width s of the wave packet was estimated 
in [13]. It was found that the anomalously high 
loss coefficient of ultracold neutrons in storage 
vessels can be explained, if the wave function of 
a neutron is represented by the de Broglie wave 
packet with the width s ~ 4 • 10~^k, where k is 
the neutron wave number. We use this estimate 
in the rest of this article. 

With the above estimate for s, if we take k^, k^ ~ 
^/u = 4 • 10^ (which is typical for thermal neu- 
trons) the angle is 7 1.3 • 10"^. This value 
is too small to measure. However, it is possible 
to design an experiment with k'^ ~ 0.1^/u, then 
7Ri 13-10"^ rad, which can be measured. 

In the case of the wave packets the change of the 

wave vector is determined by (27), Skz = 2,s"^/A;^. 
If this change is the result of the rotation by the 
angle 7, then 7 = 2s^/ k^k'^. Thus, using the same 
estimates as in the previous paragraph we obtain 
the same number for the angle 7. 

We can check this result with the help of an ex- 
periment, schematic of which is shown in fig. 3. 
In the proposed experiment, the reflection of the 
monochromatic polarized neutrons from a mag- 
netic mirror at a fixed angle is measured. After re- 
flecting from the magnetic mirror, the polarized 
neutrons go to a single crystal where the Bragg 
condition is satisfled. After reflecting from the 
crystal, the neutrons are registered by a detec- 
tor. When neutrons arc polarized parallel to the 
magnetization of the mirror, they are totally and 
specularly reflected from it. If the polarization is 
opposite, the reflection is partial and not specu- 
lar. The reflected beam from the mirror does not 
completely satisfy Bragg condition at the single 
crystal, so in order to restore the Bragg condition 
we need to rotate the crystal. Thus the experi- 
ment consists of measuring rocking curves of a 
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Figure 3. Scheme of experiment to check deviation 
of reflected neutrons from specular direction when 
reflection is partial. Monochromatic polarized neu- 
trons (pn) go through spin flipper (SF) and are re- 
flected from a magnetized mirror (MM). If SF is 
switched off then neutron polarization is parallel 
to the mirror magnetization, reflection is total and 
specular. The neutrons reflected from MM go to Si 
single crystal and after Bragg reflection to detector 
D. When SF is switched on, neutron polarization 
becomes opposite to magnetization of MM and the 
reflection is not total, thus the Bragg condition is 
not satisfied for neutrons going to Si crystal. In or- 
der to restore Bragg condition the Si crystal must 
be turned by an angle 7. The purpose of the exper- 
iment is to measure this angle. 

single crystal for two neutron polarizations. The 
centers of the two rocking curves will be shifted 
by 7. Of course, to achieve sufficient precision we 
need incident beam also monochromatized by an 
identical single crystal. More over, to avoid some 
false effects the analysis, may require reflection 
from two specially arranged single crystals. This 
can be discussed for a concrete reflectometer. 



For thermal neutrons with A = 2 A the spatial 
width of the wave packet is < l/xni, which is of 
the order of magnitude less than the G-H shift. 
However this changes, if the incident neutron in- 
stead of the wave packet is described by a beam 
of the same width. In this case we must compare 
the G-H shift to 1/scos^, which is one order of 
magnitude more than the G-H shift. 

In the case of the reflection from a layer of thick- 
ness, say 1 /im, on a substrate the G-H shift is 
in the order of 0.8 mm and is considerably larger 
even than 1/s cos 9. 

This shift is possible to measure in the experi- 
ment shown in fig. 4. If a Cadmium (Cd) plate 
restricts the area of the mirror enlightened by 
the incident beam, then the detector will detect 
only those particles, which can dive under the Cd 
plate. 

It is possible that the beam which passed under 
the Cd plate will experience multiple reflections 
from the vacuum-layer interface and from the 
layer-substrate interface. In this case the distri- 
bution of the neutrons at the position sensitive 
detector (PSD in fig. 4) will look like a diffrac- 
tion pattern. In reality, however, this pattern will 
be smeared by the non-monochromaticity and by 
scattering from the interface roughness and from 
the inhomogeneities of the film. However the con- 
tribution of scattering can be estimated and sep- 
arated, if we can measure intensity of neutrons 
scattered outside of the incidence plane. 



3.2 G-H shift at total reflection 

Let us estimate the G-H shift (15) at the to- 
tal reflection from a single interface. For thermal 
neutrons the ratio k^/ at the total reflection is 
about 400. Thus the G-H shift is 400 times larger 
than penetration depth which can be es- 

timated to be ^ 10~^ cm. Therefore, for 

thermal neutrons the G-H shift is in the order of 
8 fim. 



3. 3 G-H effect and coherence 

Reflection from two or more interfaces presents 
a problem, which is worth while to carefully in- 
vestigate. In such a reflection we flnd coherence, 
incoherence and a transition between them. 

Indeed, expression (31) is the result of a coher- 
ent superposition of plane waves instantaneously 
reflected multiple (infinite) times from two in- 
terfaces. This expression is obtained as a sum of 
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Figure 4. Multiple reflection from two interfaces 
of sufficiently thick film, when shifted wave fields 
on the mirror surface do not overlap, can be repre- 
sented as multiple incoherent reflection. In that case 
a diffraction pattern can be observed at the position 
sensitive detector PSD. Every maximum will corre- 
spond beam, which escaped the film after consecu- 
tive reflections from the substrate. 



geometric series 



J2 [e^'^^^^riae— rio 



e^^^i'^nse^^^iV, (36) 



where every term represents a possible path of 
the wave between the moment it enters the mir- 
ror and the moment it exits the mirror. The wave 
entering the mirror is defined by the vacuum- 
to-layer transmission amplitude toi — 2?^/ {qz + 
g^,), and the wave exiting the mirror is defined 
by layer- vacuum transmission amplitude tio = 
2q'^/ {qz+q'z)- If we sum these multiple reflections, 
and use the resulting reflection amplitude (31) 
to find refiected wave packet, wc obtain a single 
exit point for the packet. If wc treat every term 
of the sum separately and use series (36) to cal- 
culate reflection of the wave packet, we obtain 
many reflected wave packets each with different 
exit point. So the single G-H shift transforms into 
multiple shifts. The coordinate of the n-th exit 
point is proportional to nd. So the question is, 
what the multiple reflections are in reality: are 
they coherent sum of plane waves or a sequence 
of incoherent reflections? If they arc the coher- 
ent sum, which can be expected, because of ob- 
servable interference, then what do multiple exit 
points mean? 



It is interesting to check whether different re- 
flected wave packets interfere with each other, or 
the interference is present only when the spots 
of the consecutively reflected waves on the in- 
terface overlap (and is absent when they do not 
overlap). If they do not interfere when the spots 
do not overlap, then multiple refiections become 
incoherent. If we can measure the reflection with 
continuously changing overlapping, we can see 
how coherence transforms into incoherence. 

To check the phenomenon mentioned in the pre- 
vious paragraph, the experiment depicted in flg. 
4 can be used. If the reflected beam is originated 
at one point close to the one of incidence, then 
it will be blocked by a cadmium shutter, and 
no neutrons will be detected by the PSD. It will 
prove the coherent summation of the reflected 
amplitudes. However, if the multiple reflections 
proceed according to (36), the distribution of the 
neutrons on the PSD will be similar to a diffrac- 
tion pattern. 

There is also a different way to check the transi- 
tion of the coherence into the incoherence. Sup- 
pose the reflection from the substrate is not total, 
i.e. ri2 is a positive real number less than unity. 
In that case the full reflection amplitude (31) is 
not a unit complex number 



R{qz 



roi + ri2exp(2ig^i(i) 



1 + roiri2 e-x.^{2iqzid) 
\R{qz)\e-xjg{i^{qz)). 

Here 



(37) 



\R\ = 



rgi + 2roiri2 cos{2qzid) + rl^ , . 

-2--2- < i, (38j 



^ 1 2roiri2 cos{2qzid) r§irf2 



and ip{qz 



arcsm 



arcsm 



ri2 s\n{2qzid) 



Irli + 2roiri2 cos(2g^id) + rfa. 



roiri2 sni{2qzid) 



Y^l 2roiri2 cos(2g^i(i) + rlirl2 , 



(39) 
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Figure 5. Experimentally measured reflection curve 
for a film on a substrate [11], and its theoretical 
fitting 

If wc expand if{qz) and ln(|i?(g^)|) around the 
point qr|| = fey, then the G-H shift is a com- 
plex number. Real part of it defines the shift, 
and imaginary part defines the deviation of the 
reflected beam from the specular direction. 

A refiectometer can be used to measure the re- 
flection amplitude dependence on energy of the 
incident neutrons for a given incidence angle. 
Example of such dependence is shown in flg. 5 
(reprinted from [11]) . The most interesting fea- 
ture of this curve is the contrast of oscillations. 
Period of oscillations is related only to the width 
of the fllm, but the contrast depends on many 
factors, such as smearing of interfaces because of 
the surfaces roughness [14] , losses because of ab- 
sorption, scattering on roughnesses and inhomo- 
geneities, and possible contribution of incoherent 
reflection. 

In the experiment [11] the contribution of the 
described factors was found by fitting the theo- 
retical curve (solid line in fig. 5) to experimental 
points. The smearing (for theory see, for exam- 
ple [14]) was described by two parameters related 
to two interfaces. The contribution of the inco- 
herent reflection was described by the coherence 
length parameter. The smearing parameter char- 
acterizes the height of roughness. The coherence 
length characterizes a certain separation of en- 
ter and exit spots of the neutron wave function, 
at which the coherent and incoherent contribu- 
tions to the reflection coefficient become equal. 



The fitting in the experiment [11] was done with 
the coherence length and without it, and it was 
found that of the fitting improves consider- 
ably, if the coherence length is included. 

We show now how the coherence length, denoted 
^0, was included in theoretical description of the 
reflection curve shown in flg. 5. 

Let us introduce a function /(x), which monotonously 
increases from zero at a; = to one at a; = oo, 
and represent the reflection coefficient as the 
sum 



f{x))\Rc\' + f{x)Ri, 



(40) 



where Rc is coherent reflection amplitude, and 
Ri is the incoherent reflection coefficient. If a; = 
^x/^Oi then for small x, i.e. for a small separation 
compared to the coherence length ^0, the re- 
flection is completely coherent. For a large sep- 
aration compared to the coherence length ,^05 
the reflection is completely incoherent. In [11] the 
error function $(Cx/\/2Co) is chosen to be f{x). 
Thus the reflection coefficient is deflned as 



+ 



\R\\k^) = \Rc\\kj 



Ri{k±)<^ 



where kj_ is the component of the incident wave 
vector normal to the mirror, Rc is given by (31), 
and Rj is the incoherent reflection coefficient 




Ri 



Ri + R2 ~ 2R1R2 
1 — R1R2 



Here R12 are reflection coefficients from the two 
interfaces. In the case of ideal interfaces 



Ri 



k± — k 



i± 



k± + ki± 



R2 



ki± — k 



2± 



ki± + k2± 



In the case of smeared interfaces, when smearing 
is described by the Eckart potential [14], 



u{x) 



u 



1 + exp(x/a) 
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with smearing parameter a, 

^ _ sinh{7rai{k± - ki±)) ^ 
sinh(7rcri(/c_L + ki±)) ' 

^ _ sinh(7ro-2(fei± - /c2±)) ^ 
^ sinh(7r(T2(A;i_L + A;2±)) 
where (Ti^2 are smearing parameters of two inter- 
faces. 

The fitting of the experimental data revealed 
that the coherence parameter is equal to 1.5 
mm. Further research has to be done to under- 
stand the meaning of this coherence parameter, 
and how it relates to the width of the neutron 
wave packet. 



4 Conclusion 

We discussed the G-H effect in case of neutrons, 
however almost everything said above is applica- 
ble also to x-rays. The most important difference 
is that we do not know how to estimate the wave 
packet width for x-rays. The experiments similar 
to those proposed in this article may give at least 
an upper limit of this very interesting parameter. 



Acknowledgement 

The author is grateful to R.Gahler for dis- 
cussions, to F.Radu for useful remarks, to 
F.V.Ignatovich for help, and to E.P.Shabalin for 
support. 

References 

[1] A. Lakhtakia. Positive and Negative Goos- 
Hanchen Shift and Negative Phase- Velocity 
Mediums (alias Left-Handed Materials). 
|arXiv:physi cs/0305133 vl 30 May 2003. 

[2] F. Goos, H. Hanchen. Ein neuer und 
fundamentaler Versuch zur Totalreflexion, Ann. 
der Phys. (1947) 1 Hf. 6, S. 333-346. 



[3] F. Goos, H. Hanchen. Neumessung des 
Strahlversetzungseffekten bei Totalreflexion 
Ann. der Phys. (1949) 5 Hf. 5, S. 251-252. 

[4] I. V. Shadrivov, A. A. Zharov, Yu. S. 
Kivshar. Giant Goos-Hanchen effect at the 
reflection from left-handed metamaterials, 
arXiv: physics/030 5032 ^1 9 May 2003. 

[5] M. Maaza, B. Pardo, On the possibility to 
observe the longitudinal Goos-Hanchen shift 
with cold neutrons. Optics Communications 
(1997), 142, 84-90. 

[6] P. D. Kukharchik, V. M. Serdyuk, I. A. 
Titovitskiy, Total reflection of a Gaussian light 
beam. ZhTF 69, (1999) 74. (see Sov. Phys. 
J.Tech.Phys.) 

[7] C.G.Darwin. Phil. Mag. (1914), 27, 675-686. 

[8] G.H.Godfrey Austral. J. Phys. (1957), 10, 1-15. 

[9] L.G. Parratt, Phys. Rev. 95, 359, 1954. 

[10] V.K. Ignatovich, F.V. Ignatovitch, D.R. 
Andersen, Algebraic description of multilayer 
systems with resonances. Particles and Nuclei 
Lett. 3 [100], pp. 48-61, 2000. 

[11] D.A. Korneev, V.I. Bodnarchuk, S.P. 
Yaradaikin, V.F. Peresedov, V.K. Ignatovich, 
A. Menelle, R. Gaehler. Refelectometry study 
of the Coherent Properties of Neutron. Physica 
B 276-278, p. 973-974, 2000. 

[12] L. de Broglie, Non-Linear Wave Mechanics: 
A Causal Interpretation (Elsevier, Amsterdam, 
1960) 

[13] V.K. Ignatovich, M. Utsuro, Tentative solution 
of UCN problem Phys. Lett. A 225, 195-202 
(1997) 

[14] D. A. Korneev, V. K. Ignatovich, S. P. 
Yaradaikin, V. I. Bodnarchuk, Reflection 
of neutrons from potentials with smeared 
interfaces. JINR P4-2002, Dubna, 2002. 



12 



